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Abstract: We introduce and initiate the investigation of a general class of 4d, M = 1 



>: 

quiver gauge theories whose Lagrangian is defined by a bipartite graph on a Riemann 
surface, with or without boundaries. We refer to such class of theories as Bipartite 
Field Theories (BFTs). BFTs underlie a wide spectrum of interesting physical systems, 
including: D3-branes probing toric Calabi-Yau 3-folds, their mirror configurations of 
CN| ! D6-branes, cluster integrable systems in (0+1) dimensions and leading singularities 

in scattering amplitudes for M = 4 SYM. While our discussion is fully general, we 
focus on models that are relevant for scattering amplitudes. We investigate the BFT 
perspective on graph modifications, the emergence of Calabi-Yau manifolds (which 
arise as the master and moduli spaces of BFTs), the translation between square moves 
in the graph and Seiberg duality and the identification of dual theories by means of 
the underlying Calabi-Yaus, the phenomenon of loop reduction and the interpretation 
of the boundary operator for cells in the positive Grassmannian as higgsing in the 
BFT. We develop a technique based on generalized Kasteleyn matrices that permits 
an efficient determination of the Calabi-Yau geometries associated to arbitrary graphs. 
Our techniques allow us to go beyond the planar limit by both increasing the number 
of boundaries of the graphs and the genus of the underlying Riemann surface. Our 
investigation suggests a central role for Calabi-Yau manifolds in the context of leading 
singularities, whose full scope is yet to be uncovered. 
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1. Introduction 

In this paper we introduce and initiate the investigation of a general class of 4d, M = 1 
quiver gauge theories whose Lagrangian is defined by a bipartite graph on a Riemann 
surface, with or without boundaries. We refer to such class of theories as Bipartite 
Field Theories (BFTs). 

The motivation for studying these theories follows from the fact that they underlie 
a wide range of interesting physical systems. 1 One of their first appearances has been 
in the context of D3-branes probing toric Calabi-Yau 3-folds in Type IIB string theory. 
The superconformal field theory on the worldvolume of the D3-branes is indeed encoded 
by a bipartite graph on a 2-torus [0. When investigating the Type IIA configurations 
of D6-branes related to the previous setup by mirror symmetry, bipartite graphs on 
higher genus Riemann surfaces surprisingly emerge as fundamental objects [[|. More 
recently, bipartite graphs on a 2-torus have been shown to give rise to an infinite class 
of quantum mechanical integrable systems [[J. 

In the latest addition to this string of applications, bipartite graphs have been 
used in the context of scattering amplitudes in quantum field theory (QFT) [|J. In 
recent years, we have witnessed tremendous progress in our understanding of scattering 
amplitudes in gauge theory, most notably for M = 4 super Yang-Mills in the planar 
limit. These developments were originally triggered by Witten's twistor string || and 
have resulted in efficient tools for the computation of scattering amplitudes at tree 
level, such as CSW diagrams || and BCFW recursion relations [[7| |[ |9|, |10| , and loop 
level [[IT], |12|, 13, 14]. At the same time, the hidden dual superconformal symmetry 



of planar M = 4 SYM was unveiled [O, [161, and it was realized that superconformal 



and dual superconformal symmetries combine to give rise to the infinite dimensional 



1 Here we do not mention the obvious, and extremely interesting, applications of bipartite graphs 
to condensed matter physics. 
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Yangian symmetry fL7 |. In order to better understand the role of this symmetry, it 
was suggested one should focus on the leading singularities of scattering amplitudes, 
which in turn arise as residues of a contour integral over the Grassmannian |L8). At 
a fundamental level, the new insights have led to the idea that a new formulation of 
QFT might exist, displaying its otherwise hidden simplicity by abandoning manifest 
locality and unitarity in favor of making the infinite Yangian symmetry explicit. The 
fundamental structure behind this new formulation might take several tightly related 
disguises: the Grassmannian, algebraic geometry, certain graphs or, as we advocate in 
this paper, certain quiver gauge theories. 

For all the systems mentioned above, BFTs are not merely a different interpretation 
of the same underlying graph. In fact, every statement and computation in any of these 
systems has a counterpart in the corresponding BFT. As usual, having an alternative 
perspective on the same physics, in this case in terms of a gauge theory, is extremely 
valuable. It not only allows us to understand known facts in a new light, but it provides 
intuition, suggests new approaches and might eventually become useful for answering 
new questions. The dynamics, duality and natural connection to geometry in the form 
of moduli spaces of the gauge theory have fruitful applications to the other systems. 

The appearance of quivers in connection to all these systems is not a coincidence. 
To some extent, it can be understood as the result of a powerful mathematical structure 



that underlies all of them: cluster algebras |19| . Quivers indeed provide the natural 
physical arena for cluster algebras. 

In this paper we will study BFTs in full generality, but our examples will mainly 
focus on graphs with boundaries, which is the sub-class that has been less explored 
so far and the one that is relevant for leading singularities in scattering amplitudes. 
Our results will also shed light in the relatively unexplored area of the combinatorics 
of bipartite graphs on general Riemann surfaces with boundaries. 

This paper is organized as follows. BFTs are introduced in Section Section 
[3] explains the BFT interpretation of various modifications of the underlying graph. 
Section f| provides an overview of several areas in which BFTs arise. In Section |5|, 
we generalize the method based on the Kasteleyn matrix for the determination of 
perfect matchings to deal with bipartite graphs with boundaries. Section [| discusses 
two toric Calabi-Yau (CY) manifolds that are associated to any BFT, its master and 
moduli spaces. Section ^ collects explicit examples of BFTs. The examples presented 
are related to leading singularities and go beyond the planar limit in two directions, 
increasing both the number of boundaries of the graph and the genus of the underlying 
Riemann surface. Section |8] studies the action of Seiberg duality on BFTs, which 
translates to square moves in the underlying graph. We analyze both planar and non- 
planar graphs, explicitly illustrating that the moduli is invariant under Seiberg duality 
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and can be efficiently exploited for identifying theories related by square moves. In 
Section |9|, we explain how the number of loops in certain diagrams can be reduced and 
explain this process in terms of the dynamics of the associated BFT. The equivalence 
between different multi-loop diagrams has a striking manifestation in terms of a single 



underlying Calabi-Yau manifold, arising as the moduli space of the BFTs. Section |373 
explains how the boundary operator on a cell in the positive Grassmannian maps to 
the Higgs mechanism in the corresponding BFT. We show how a new bipartite graph, 
related to the original one by the untwisting map on zig-zag paths, is an efficient tool 
for identifying consistent higgsings. We conclude and summarize some open questions 



for further research in Section 11 



Note added: while this paper was being finalized, we became aware of [p7| , which has 
some overlap with this work. 



2. Bipartite Field Theories 

In this section we introduce the concept of a Bipartite Field Theory. A BFT is a 4d, 
M = 1 quiver gauge theory whose Lagrangian is defined by a bipartite graph G on a 
Riemann surface E, which can contain boundaries. BFTs are natural generalizations 
of toric quivers, which are defined by bipartite graphs without boundaries on T 2 Q. 




Figure 1: An artistic representation of a bipartite graph on a Riemann surface defining a 
BFT. 



We can separate the nodes of G into internal and external (those on boundaries) 
ones. G is bipartite if nodes can be colored in black and white such that: 
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1) Every edge in G connects nodes of different colors. 

2) Every boundary node is connected to a single edge. 

2.1 The Dictionary 

Let us now describe the translation between the elements of the bipartite graph on E 
and a 4d, Af = 1 quiver gauge theory. 

• Faces : to every face in the graph, we associate an SU(N) group. 2 We can identify 
faces as internal or external. Internal faces are those whose entire perimeter is 
given by edges in G. On the other hand, external faces are those in which part 
of its perimeter overlaps with the boundary of the graph. We show examples of 
both types of faces in Figure £| 

• Edges: every edge is identified with a chiral multiplet transforming in the bi- 
fundamental representation of the two faces it separates. The bipartiteness of 
G introduces a natural orientation of the bifundamentals dual to edges which, 
without loss of generality, can be taken to be oriented clockwise around white 
nodes and counterclockwise around black nodes. 

• Nodes : every internal node corresponds to a monomial in the superpotential, 
given by the product of all chiral fields associated with the edges terminating on 
it. The valence of the node, i.e. the number of edges terminating on the node, 
corresponds to the order of the superpotential term. We assign signs to superpo- 
tential terms such that white and black nodes correspond to plus and minus signs, 
respectively. External nodes do not have any superpotential interpretation, and 
simply follow from the existence of edges that have one endpoint on an internal 
node and terminate on the boundary. 

Due to the bipartiteness of G, the number of edges around an internal face is even. 
Furthermore, the numbers of white and black nodes around an internal face are equal. 
In contrast, the number of edges on the perimeter of an external face can be either 
even or odd. Given the interpretation of edges as bifundamental chiral fields with an 
orientation dictated by the color of nodes we conclude that, from a quiver perspective, 

2 In fact, each face corresponds to a U(N) group. When gauging some of them, the U(l) pieces are 
IR free and become global symmetries at low energies. 
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internal faces have an equal number of incoming and outgoing bifundamentals arrows, 
which in turn implies that they are free of anomalies. As a result, the SU(N) groups 
associated to internal faces can be consistently gauged. From now on, we include in 
the definition of a BFT the fact that internal faces are gauged, while external faces 
correspond to global symmetry groups. 3 ' 4 



This interpretation of faces is rather natural in light of similar systems that admit 
an interpretation as brane configurations in string theory. In such systems, faces corre- 
spond to stacks of D-branes suspended from a web of branes associated with the graph 
G, which extends over dimensions transverse to the ones in which the gauge theory 
lives. Internal faces have a finite extension along the directions transverse to the field 
theory ones and hence give rise to gauge symmetries. External faces can instead have 
an infinite extension in these directions and thus lead to global symmetries. BFTs 
associated to graphs with no boundaries on T 2 indeed arise as configurations of D5 and 
NS5-branes The question of whether general BFTs can arise as systems of branes 
is a very interesting one, but we postpone it for future studies. 

According to the map between faces and gauge or global symmetry groups, edges 
can correspond to: bifundamentals of the gauge group (when they sit between two 
internal faces), fundamental or antifundamental flavors (internal/external) or gauge 
singlets transforming in a bifundamental representation of the global symmetry group 
(external/external). 

3 Following the previous discussion, there exists an alternative natural gauging, in which a basis for 
all closed loops in the graph whose perimeter is made out of internal edges is gauged. This alternative 
becomes relevant for non-planar graphs and is useful in some contexts such as leading singularities. 
While we will not consider this gauging any further, the resulting theories can be studied with exactly 
the same tools we develop in this paper. 

4 Internal faces might contain punctures inside them. A possible way of distinguishing such faces 
at the level of the BFT is by not gauging the corresponding anomaly-free SU(N) group. While in 
this paper we will gauge the symmetries associated to all internal faces, it is important to keep this 
possibility in mind. 




nternal 



External 



Figure 2: Examples of internal and external faces in a BFT. 
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Figure 3: A bipartite graph and its quiver dual. Blue nodes in the quiver represent gauge 
groups and green squares represent global symmetry groups. Plaquettes in the quiver corre- 
spond to superpotential terms. 



The correspondence we have introduced implies that every BFT has a quiver dia- 
gram living on E, which is indeed dual to the bipartite graph G as illustrated with an 
example in Figure |3|. Faces, edges and nodes in G are mapped to symmetry group nodes 
(gauged or global), arrows and plaquettes (superpotential terms with sign determined 
by their clockwise or counterclockwise orientation) in the quiver, respectively. Below 
we summarize the dictionary between graphs and BFTs. 



Graph 


BFT 


Internal face (2n-sided) 


Gauge group with n flavors 


External face 


Global symmetry group 


Edge between two faces i and j 


Chiral multiplet in the bifundamental 
representation of the groups i and j. The 
orientation of the corresponding arrow is such 
that it goes clockwise around white nodes and 
counterclockwise around black nodes 


A;-valent node 


Monomial in the superpotential involving k chiral 
multiplets. The signs of the terms are (+/-) for 
(white/black) nodes 



Let us briefly comment on the scale dependence properties of BFTs. A natural 
representation of bipartite graphs is given by the so called isoradial embedding. We refer 
the reader to |20) for details about this construction. Interestingly, in this embedding 
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it is possible to map the scaling dimensions of chiral fields to the angles subtended by 
the corresponding edges. Vanishing of the beta functions for gauge and superpotential 



couplings then translates into local flatness of the graph jy, |2TJ. As a result, BFTs on 
curved Riemann surfaces are not conformal. 

In the remainder of this section, we discuss two important concepts in the study 
of bipartite graphs: perfect matchings and zig-zag paths. 



2.2 Perfect Matchings 

An almost perfect matching p is a subset of the edges in G such that: 

• Every internal node is the endpoint of exactly on edge in p. 

• Every external node belongs to either one or zero edges in p. 

Almost perfect matchings can be regarded as regular perfect matchings of a larger 
bipartite graph that contains G and extends it beyond its boundaries. For brevity, we 
simply refer to them as perfect matchings in what follows. Perfect matchings connect 
bipartite graphs to gauge theory, toric geometry and integrable systems, as we explain 
later. 

The map between chiral fields in the quiver X i: equivalently edges in G, and perfect 
matchings p^ is given by 

c 

where c is the total number of perfect matchings, and P^ is equal to 1 if the edge in the 
bipartite graph associated to the chiral field Xi is contained in p^ and zero otherwise 
i, i-e. 

flifJQe^ 

In Section |5|, we will introduce efficient methods for determining the matrix P. It will 
play an important role in Section 0, when computing the moduli spaces of the BFTs. 5 



5 Throughout the paper, we will also use a notation for chiral fields involving two subindices, ex- 
plicitly indicating the gauge groups under which they are charged, as opposed to the single subindex 
notation in (2.2). We arc confident that these two alternative notations will not generate any confusion. 
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2.3 Zig-Zag Paths 

Zig-zag paths, also denoted alternating strands, are oriented paths in a bipartite graph 
that alternate between turning maximally right and maximally left at every node. They 
can be efficiently implemented in terms of a double line notation for edges [[| , in which 
two zig-zag paths go over every edge in opposite directions, crossing at the middle 
point. Figure [| shows the zig-zag paths for the example in Figure 




Figure 4: Zig-zag paths in double line notation. This example has five zig-zag paths going 
through external legs and an internal one, shown in purple. 



Zig-zag paths play a prominent role in the study of bipartite graphs for various 
reasons. First, it is possible to reconstruct the graph G from knowledge of its zig-zags. 
There is an edge for every intersection between a pair of them, and black (white) nodes 
correspond to disks with clockwise (counterclockwise) oriented boundaries. Relative 
motion of the zig-zags results in different graphs, which translate to dual BFTs as we 
will briefly discuss in Section^ (see also In the case of graphs without boundaries 

on T 2 , zig-zag paths are directly connected to the Calabi-Yau geometry of the moduli 



space of the associated BFT |fj, p2| . They are the legs in the (p, q)-weh [^3], ^4], |25 
dual to the corresponding toric diagram. Whether there also exist a direct link between 
zig-zags paths and moduli spaces in the generalized context of arbitrary BFTs is an 
extremely interesting question that we plan to revisit in the future. 

Consistency: we will restrict to graphs that do not contain self-intersecting zig-zag 
paths. The existence of such paths has, in some sub-classes of BFTs, been linked to 
inconsistencies of the field theory [^2j . For brevity, we will refer to the resulting graphs 
and BFTs as consistent. A full study of the consequences of self-intersecting zig-zags 
on general BFTs is beyond the scope of this paper. We will often explicitly determine 
the zig-zags of graphs, showing this condition is met. 
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Untwisting Map: the action of the untwisting map is schematically shown in Fig- 
ure |5[ It interchanges: 



G on £ G on E 

zig-zag path <H- face 
face <r¥ zig-zag path 

The graph G on £ resulting from untwisting can be interpreted as a new bipartite 
field theory, which we denote BFT. The concept of BFT will be used in Section |10| for 
keeping track of zig-zag paths under certain modifications of the graph. 



A 2 Aj 



Bj B 2 
Figure 5: The untwisting map. 



Deligne Permutations: for bipartite graphs with boundaries, it is possible to intro- 
duce the concept of Deligne permutations jp, which correspond to a bijection fd{bi) = 
bj that maps every boundary point 6, to another boundary point bj. 6 More con- 
cretely, Deligne permutations are in one-to-one correspondence with zig-zag paths going 
through external legs as follows 

fd(bj) = bj -<=>- The zig-zag path starting at 6j ends at bj (2.3) 
For example, for the model in Figure |] we have 

f d (l) = 4, / d (2) = 5, / d (3) = l, 

f d {A) = 2, / d (5) = 3. 1 • ] 



3. A BFT Perspective on Graph Modifications 

In this section we discuss the BFT interpretation of various possible modifications of 
the bipartite graph. 

6 Othcr works in which bipartite graphs have been discussed in conjunction with certain permuta- 
tions and Calabi-Yau manifolds include [^9[ [30| . 
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3.1 Reduction to 2 and 3-valent Graphs 

Nodes with valence k in the graph correspond to order k superpotential terms in the 
BFT. In particular, 2-valent nodes correspond to mass terms. Due to the specific struc- 
ture of the BFT superpotentials, integrating out massive fields has a simple graphical 
implementation: it corresponds to condensing the two nodes at both endpoints of the 
2-valent one, as illustrated in Figure |. 




Pi(Xi) P 2 (Xj) Pi(Xi) x P 2 (Xi) 

Figure 6: Integrating out massive fields associated to 2-valent nodes corresponds to con- 
densing nodes on the graph. 



The origin of the condensation can be understood as follows. The superpotential 
takes the form 

W = X 1 P 1 (X i ) + X 2 P 2 (X t ) - X,X 2 + . . . (3.1) 

where we have identified all terms in the superpotential containing X± and X 2 . P\(Xj) 
and P 2 (Xi) are products of bifundamentals fields. The equations of motion for the 
massive fields are 

d Xl W = 0^X 2 = P 1 (X i ) and d X2 W = & X 1 = P 2 (JQ). (3.2) 
Plugging them back into (|3.1| ), the terms involving X\ and X 2 are replaced by 



W = P l (X t )P 2 (X t ) + ... J (3.3) 

which is precisely the interaction associated to merging the nodes at both sides of the 
mass term. I.e. two terms of order k\ and k 2 are combined into a single one of order 
h + k 2 - 2. 

Inverting this process, we can reduce the order of superpotential terms by inserting 
2-valent nodes, i.e. by integrating-in massive fields. Iterating this process, it is possible 
to reduce any fc-valent node to (k — 2) 3-valent and (k — 3) 2-valent ones. This process is 
clearly not unique, although the low energy physics is independent of how we perform 
it. Figure [7| shows a possible decomposition of a 5-valent node. 
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Figure 7: A possible decomposition of a 5-valent into three 3-valent ones by the insertion of 
two 2-valent nodes. The field theory counterpart of this process corresponds to integrating-in 
massive fields. 



Prom Bicolored to Bipartite Graphs 

Bicolored graphs that are not bipartite, i.e. graphs that contain edges connecting nodes 
of the same color, can be studied with the same tools discussed in this paper. Plabic 
(i.e. planar bicolored) graphs are examples of this class of models [p6 . 



Whenever we encounter a graph with an edge connecting two nodes of the same 
color, we will interpret it as a bipartite one by introducing a 2-valent node of the 
opposite color in the middle of this edge. The corresponding massive fields can then be 
integrated out, resulting in the merging of the two original nodes. This field theoretic 
interpretation leads precisely to the two equivalent procedures that appear in the math 
literature for turning bicolored graphs into bipartite ones: introducing an intermediate 
2-valent node or condensing the two nodes of the same color. 
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Figure 8: A bicolored graph containing an edge between nodes of the same color can be 
turned into a bipartite one by introducing an intermediate 2-valent node or by node conden- 
sation. 



3.2 Seiberg Duality 

Seiberg duality |27j is a remarkable property of Af = 1 gauge theories that implies 
a full equivalence between two different theories (typically referred to as electric and 
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magnetic) in the low energy limit. It plays an essential role in our understanding of 



the low energy dynamics of SUSY gauge theories [28|. Being such a central concept in 



SUSY field theories, Seiberg duality has, in the BFT context, a natural implementation 
in terms of the underlying bipartite graph. 

We will focus on Seiberg dualities in which the electric and magnetic theories are 
BFTs on the same Riemann surface. This constraint implies that we can only dualize 
Nf = 2N C gauge groups, where Nf = N, the common rank of all gauge and global 
symmetry groups. This class of gauge groups corresponds to internal square faces 
in the graph. Let us briefly discuss the field theory side of the duality, focusing on 
the dualized gauge group, while considering all others as expectators. The electric 
theory has SU(N C ) gauge groups and flavors Qi, Q 2 in its fundamental representation 
and Qi and Q 2 in the antifundamental representation. More precisely, these flavors 
are in fact bifundamentals fields in the BFT but, for brevity, we omit discussion of 
the transformation properties under the additional symmetry groups. The dual gauge 
group is SU(N C ), where the dual rank is N c = Nf — N c = N c , i.e. it is equal to the 
original one. In addition, the duality implies the following transformations: 



• Replace the electric quarks Q%, Q 2 , Q\ and Q 2 by magnetic quarks q%, q 2 , q± and 
q 2 transforming in the conjugated (bifundamental) representations. 

• Introduce meson fields My, which are composite from the perspective of the 
electric theory, i.e. My = QiQj. Meson fields are singlets of the magnetic gauge 
group and transform in bifundamentals representations of the other groups, which 
are either gauged or global. 

• Introduce cubic superpotential couplings between the dual quarks and the mesons 

AW = ^2q i M ij q j , 

ij 

and re-express any product of the electric quarks in the original superpotential 
in terms of the mesons. 



In Figure |9| we present the simple graph transformation that implements Seiberg 
duality. It was originally discovered in [0 in the context of BFTs without boundaries 
on T 2 but extends to generic BFTs without changes. It is easy to see that this modi- 
fication of the graph automatically implements the three points discussed above. It is 
often referred to as urban renewal, spider move or square move. If 2-valent nodes are 
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generated during this process, the corresponding massive fields can be integrated out 
as discussed in Section |37T| . Repeating this operation twice on the same face obviously 
returns to the original graph. 



0' 



M 21 



Figure 9: Seiberg duality maps to square moves in the graph. 



The effect of Seiberg duality on zig-zag paths is shown in Figure [10]. It corresponds 
to a reorganization of the four zig-zag paths passing through the dualized square, in 
which zig-zags associated to opposite corners are pairwise interchanged. 




■=> 




Figure 10: The action of Seiberg duality on zig-zag paths. 



Seiberg duality is a full equivalence between the electric and magnetic theories in 
the IR limit. This implies agreement in their behavior under relevant deformations, 
matching of their moduli spaces, etc. The field theoretic perspective on bipartite graphs 
thus suggests natural invariants under square moves. This intuition can be exploited in 
the context of any other interpretation of the same graphs, such as leading singularities 
in scattering amplitudes. We will use this approach in Sections |] and 



3.3 Higgsing 

Edge removal is another natural operation on graphs. From a BFT viewpoint, it 
translates to the scalar in the corresponding bifundamentals chiral multiplet acquiring 
a non-zero vacuum expectation value (vev) jjj]. The two faces on each side of the 
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removed edge are combined into a single one as shown in Figure [11]. Depending on the 
type of faces that are merged, we can have three different situations, which have the 
following BFT interpretation: 



• Internal-internal: higgsing of the corresponding SU(N) x SU(N) piece of the 
gauge group down to the diagonal SU (N) subgroup. 

• Internal-external: color-flavor locking of the SU(N) gauge x SU(N)gi b a i sym- 
metry associated to the faces. 

• External-external: spontaneous breaking of an SU(N) x SU (N) subgroup of 
the global symmetry to the diagonal subgroup. This process results in massless 
Goldstone bosons. 



If 2-valent nodes are generated in this process, the corresponding massive fields can be 
integrated out. Throughout the paper we will be interested in preserving external legs, 
which map to scattered particles, so we will no longer consider the third option. 



cr^v> 



1/2 



Figure 11: Removing an edge in the graph corresponds to turning on a vev for a bifunda- 
mentals scalar, resulting in the merging of two faces. 



In Section ID we will investigate the connection between higgsing and the boundary 
operator in leading singularities and introduce an implementation of it that efficiently 
keeps track of zig-zag paths, ensuring that we only perform higgsings that result in 
consistent graphs. 



4. BFTs Everywhere 

BFTs and their dynamics play an important role in a wide spectrum of interesting 
physical systems. In this section we present a brief overview of some of them. 
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4.1 D3-Branes over Toric Calabi-Yau 3-Folds 



The 4d, M = 1 superconformal field theory (SCFT) arising on a stack of D3-branes 
probing a CY 3- fold is a BFT on a 2-torus without boundaries The two fundamental 
directions of the T 2 correspond to a U(l) 2 flavor symmetry that follows from isometries 
of the toric CY. The remaining U(l) isometry translates to the R-charge of the gauge 



theory. Figure [12] shows an example, corresponding to the BFT on D3-branes at the 
complex cone over F [0. 7 




Figure 12: Bipartite graph associated to the 4d, N = 1 SCFT on D3-branes over the 
complex cone over Fq . Opposite sides of the square are identified to form a T 2 . 



The correspondence between this class of SCFTs and bipartite graphs has been 
extremely fruitful, fully answering the question of which gauge theory is associated to 
a given toric Calabi-Yau and vice versa. 

4.2 Mirror Symmetry 

Consider the configuration discussed in the previous section, with D3-branes on a toric 
singularity with characteristic polynomial P(zi,Z2) = J2 a ni,n 2 z T z 2 2 i where (nj,^) 
runs over points in the toric diagram. On the D3-branes, we obtain a BFT theory 
described by a bipartite graph G on a 2-torus. 

The mirror manifold is given by P{z\,z-i) = W, W = uv. Let us consider the 
Riemann surface S sitting at W = 0. The genus and number of punctures of E 
are given by the number of internal points and the perimeter of the toric diagram, 
respectively. The original configuration of D3-branes is mapped to a set of intersecting 
D6-branes in the mirror, with one type of D6-branes for each gauge group in the quiver 
and bifundamentals chiral multiplets arising at their intersections f2|. The non-trivial 
information about the mirror configuration is encoded in a new bipartite graph G on 

7 In fact there is another BFT associated to D3-bran.es on the same geometry, which is connected 
to this one by Seiberg duality as discussed in Section g. 
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E, in which the D6-branes arise as zig-zag paths. The graph G is obtained from G 
by applying the untwisting map. We can associate a new BFT gauge theory to G, as 
mentioned in Section |2.3| . The BFT has been referred to as the specular dual in 



Figure |13] shows G obtained by untwisting the Fq model in Figure |12[ In this case, £ 
is a 2-torus with four punctures. 





Figure 13: The tiling of the mirror obtained by acting with the untwisting map on the Fq 



model in Figure 12. The resulting Riemann surface on which the graph lives is a 2-torus with 
four punctures. 



The discussion in this section makes it clear that mirror symmetry naturally gives 
rise to bipartite graphs, i.e. BFTs, on Riemann surface of arbitrarily high genus, 
providing an important motivation for their study. 

4.3 Cluster Integrable Systems 

Bipartite graphs without boundaries on a 2-torus are also in one-to-one correspondence 
with an infinite class integrable systems in (0+1) dimensions, denoted cluster integrable 
systems || . The Poisson manifold of the integrable system is parametrized by loops on 
the graph, a useful basis for which is provided by the loops around faces and along the 
two fundamental directions of the T 2 . The Poisson brackets between these variables are 
dictated by the number of edges over which the corresponding loops overlap, counted 
with orientation. Different patches of the Poisson manifold, typically described by 
different graphs, are connected by cluster transformations. In [EJ, it was shown that 
the Hamiltonian and Casimir operators of the integrable system correspond to internal 
points and ratios of external points in the associated toric diagrams, respectively. We 
refer the reader to || for a detailed explanation of how the integrals of motion are 
constructed in terms perfect matchings. Figure [14| shows a graph which, by means 
of this correspondence, is mapped to the n-particle, relativistic, periodic Toda chain. 
Physical implications of this correspondence and connections to other setups realizing 



the same integrable systems have been studied in [32, 33, 34, 35 . 
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Figure 14: Bipartite graph on a T 2 associated to the n-particle, relativistic, periodic Toda 
chain. 



Motivated by generic BFTs, it is natural to ask whether bipartite graphs on Rie- 
mann surfaces other than T 2 also give rise to integrable systems. We expect it might 
be possible to extend the proof in [[J to at least graphs without external nodes that 
can be embedded in a T 2 , i.e. disks and cylinders. 



4.4 Leading Singularities 

The connection between bipartite graphs, and consequently BFTs, and the positive 
Grassmannian and leading singularities in planar M = 4 SYM is relatively more recent 
and less known to the general physics audience. For this reason, our discussion of these 
topics will be slightly lengthier. The material presented in this subsection is a brief 



summary, tailored for our specific needs, of the ideas in P, 26] and references therein, 
to which we refer the interested reader for further details. 

Leading singularities can be found by an iterative procedure. Starting from a 
loop amplitude, we compute the discontinuity at branch cuts, next we determine the 
discontinuity at branch cuts of this result, and repeat the process until we are left with 
a rational function of the kinematical invariants. For brevity, this discontinuity across 
the leading singularity is often referred to as the leading singularity. We can regard 
leading singularities as calculable well-defined data associated to gauge theories. In 
addition it has been conjectured that, for maximally supersymmetric theories, leading 
singularities are sufficient for determining the perturbative S Matrix Evidence 



supporting this proposal for Af = 4 SYM has been given in [|6|, |37|, |38|, |39 



The Grassmannian G(k, n) is the space of A>dimensional planes in n dimensions. 
Points in G(k,n) can thus be parametrized by a general k x n matrix C, whose rows 
correspond to n-dimensional vectors spanning a plane. We can take any linear combi- 
nation of the rows without affecting the plane, so we conclude that this parametrization 
has a GL(k) redundancy or 'gauge symmetry'. The positive part of the Grassmannian 
G>o(k,n) corresponds to the subspace in the Grassmannian in which the determinant 
of all A; x A; minors of C are greater or equal to zero. In what follows, n corresponds to 
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the total number of scattered particles, with k being the number of negative helicity 
ones. Amplitudes with k = 2 are known as maximally helicity violating (MHV), while 
for k > 2 they are denoted N fc_2 MHV. 

It has been proposed in Jl8[ that all leading singularities in planar M = 4 SYM 
arise as residues of the following contour integral over the Grassmannian 



(12-t)(23-(t + l))-( B l-(-l)) n^(W (4.1) 

Here . . ik) indicates the determinant of the k x k matrix made out of the i\, . . . , %k 
columns of C, i.e. the denominator consists of the determinants of the n sequential 
minors in C . Wj = (Xj, fij, fjj), j = 1, . . . , n, are the kinematic variables of the scattered 



particles in twistor space. We refer the reader to [|18] for a detailed explanation of this 
representation and of how to identify the contour integration, which determines the 
resulting leading singularity. Leading singularities correspond to certain subspaces, 
also denoted cells, of the Grassmmanian parametrized by a constrained matrix C. 

On a parallel line of development, subspaces of the Grassmannian have been shown 
to be in one-to-one correspondence with bipartite graphs P5| . More concretely, cells 
in G(k, n) are associated to bipartite graphs on a disk with n boundary points, k of 
which are black. 8 We thus have a connection between the following objects: 

Leading singularities <^ Cells in the Grassmannian Bipartite graphs 

Let us explain how to go from a bipartite graph to a cell in the Grassmannian. The 
first step is to define certain 'momentum flows' along edges of the graph, also called 
perfect orientations, which are in one-to-one correspondence with perfect matchings 
p6| . These flows are such that there are two outgoing and one incoming arrows at each 



white node and two incoming and one outgoing arrows at each black node, as shown 
in Figure [T5|. Flows go through 2-valent nodes without changing direction. Since every 



bipartite graph can be reduced to 2 and 3-valents nodes as explained in Section [3T 
these rules are sufficient for determining perfect orientations. 

The bijection between perfect matchings and perfect orientations works as follows. 
Given an edge contained in a perfect matching, we identify it with the incoming and 
outgoing arrows of the white and black nodes at its endpoints, respectively, as shown 
in Figure [16]. In Figure IT7| we present an example of a perfect matching and its corre- 



sponding perfect orientation. 

8 More generally, cells in the Grassmannian can be parametrized by plabic graphs, but it is straight- 



forward to turn them into bipartite ones, as explained in Section 3.1 
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Figure 15: Flow pattern on 3-valent nodes. 





Figure 16: There is a one-to-one correspondence between perfect matchings and perfect 
orientations. Here we indicate in red an edge in a perfect matching. 



2 3 



2 3 




8 

7 6 7 6 

(a) (b) 
Figure 17: a) An example of a perfect matching and b) its corresponding perfect orientation. 



We are ready for determining the restricted structure of C associated to a graph, 
given a perfect orientation. C is a k x n matrix in which now rows correspond to 
black external nodes and columns correspond to all external ones, i.e. negative and 
positive helicity particles are mapped to black and white external nodes, respectively. 9 
Entries for which the row and column correspond to the same node are set to 1. Entries 
associated to a pair of nodes that are not connected by an oriented path in the perfect 
orientation are set to 0. Finally, all other entries are determined in terms of edge 



weights via the so called boundary measurement [26], 



9 Thc explicit form of C depends on how we arrange nodes within the rows and columns. This 
freedom can be taken care of by the existing GL{k) 'gauge' symmetry. 
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e n*. 

':6i->6i e£P 



(4.2) 



i.e. we sum over all directed paths P starting from 6j and terminating at bj in the 
graph with a perfect orientation, and the product is over all edges e in P. When P 
has self-intersections, we have to weigh the corresponding contribution by (— i) wmd ( p ) ) 
where the winding wind(P) is the signed number of full 360° turns P makes. The 
edge weights x e are in one-to-one correspondence with the expectation values of the 
corresponding scalars in the BFT. The precise map between them will be clarified in 
Section [U], when we discuss higgsing. 

The procedure we have just outlined clearly depends on a choice of perfect orien- 
tation, equivalently on a choice of perfect matching. Different choices are physically 
equivalent. The result of this prescription coincides with the kinematical analysis of 
leading singularities, as the one presented in [EDI. For example, for the configuration 



in Figure [T7|, the C matrix becomes 



C 





1 


2 
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5 
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7 
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2 





1 





C24 





C26 


C27 





3 
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C34 





C36 


C37 





5 











C54 
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C56 


C57 
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C84 





C86 


C87 


1 



(4.3) 



The subspace parametrized by the constrained matrix C associated to a bipartite graph 
when edge weights are restricted to be R > 0, is a cell in the positive Grassmannian. 

The previous discussion makes the connection between bipartite graphs and cells 
in the Grassmannian relatively natural. White and black nodes can be interpreted as 
MHV and MHV 3-point amplitudes and, intuitively, the graph provides a picture of 
a scattering process in which external nodes represent scattered particles and internal 
faces correspond to loops. 



5. Kasteleyn Technology for General BFTs 

In this section we introduce an efficient method for finding the perfect matchings of a 
general bipartite graph, generalizing the approach based on the Kasteleyn matrix to 
graphs that might contain boundary nodes. These techniques will play an essential role 
in the efficient computation of moduli spaces. 

We begin by defining the master Kasteleyn matrix Kq, as the adjacency matrix of 
the graph in which rows are indexed by white nodes and columns are indexed by black 
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nodes, i.e. for every edge in the bipartite graph between nodes w M and h u , we introduce 
a contribution to the K 0>flv entry. We separate white nodes into two sets Wi and W e , 
corresponding to internal and external (i.e. boundary) nodes, respectively. Similarly, 
we split black nodes into Bi and B e . This separation is independent of the number 
of boundary components and of how external nodes are distributed among them. The 
individual numbers of internal and external nodes need not be the same for different 
colors. Furthermore, the total numbers of white and black nodes need not be equal, 
either. Kq takes the general form 





f 






\ 






* 


* 






* 





7 



(5.1) 



2 




Figure 18: A BFT on a disk, a) The corresponding graph has one internal and four external 
faces. It also has four white nodes (two internal and two external) and four black nodes (two 
internal and two external), b) The associated quiver diagram. 



Let us illustrate these ideas with the simple BFT shown in Figure |18|a, which is 
related to a leading singularity in the scattering of 2 negative helicity and 2 positive 



helicity gluons at 1-loop. Figure [L8|.b shows the corresponding quiver diagram, for 
which the superpotential is 



W — X 15 X 52 X 21 — Xi 3 X 32 X 2 \ + X13X34X41 — X15X54X41, 



(5.2) 



where color indices and their contractions are implicit. The master Kasteleyn matrix 
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IS 



/ 
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7 8 
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-^15 




-^52 
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X41 







X34 


3 













4 





X32 









(5.3) 



For any subsets W e ^ e i Q W e and B e ^ e i Q B e of the boundary nodes, we define the 
reduced Kasteleyn matrix: 



K^ w d ; B ^ ^ = matrix resulting from starting from K and deleting the rows 

in W et dei and the columns in B e ^ e i (5.4) 



All perfect matchings in the graph are then encoded in the polynomial 



W e del,B e del 



(5.5) 



where the sum runs over all possible subsets W e ^ e i and B e dd of the external nodes 
such that the resulting reduced Kasteleyn matrices are square. Every term in this 
polynomial, which we denote V^, is interpreted as the product of edges in a perfect 
matching. 

Let us explain in more detail the reason for the sum over reduced Kasteleyn ma- 
trices in (WM- The determinant of each K/„ r . . a _ _\ in (|5.5| ) generates all the perfect 



y Wf, t del,B ei del J 

matchings containing the edges connected to the external nodes in (W e — W e ,tfez) and 
(B e — B e ^ei)- Once again, let us show how this works in the example in Figure [18|, for 
which K is given in (|5.3|). Let us first consider W e> d e i = {3} and B etdei = {7}, which 
results in 



K, 
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X\§ 


X21 





{3},B e , de! ={7}) 
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X41 
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X32 






\ 



(5.6) 



Taking its determinant, we obtain 

detK, 



w e . 



{3},B e 



={7}) 



-X15X32X34. 



(5.7) 



This is the only perfect matching that contains the edges connected to the surviving ex- 
ternal nodes 3 and 8, and we show it in Figure |TJ^ Generically, each reduced Kasteleyn 
matrix can give rise to multiple perfect matchings. 
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Figure 19: Perfect matching generated by detK^ w ^ ;= { 3 } b d ;={7})' Edges m the perfect 
matching are indicated in red. 



Computing the full polynomial V in (|5.5|) , we obtain 

V = ^13^15 — ^21^41+^32^41^52—^15^32^34+^21^34^54—^13^52^54+^32^34^52^54, 

(5.8) 



which corresponds to the seven perfect matchings in Figure [19]. Given the definition 
in (|2.2|) , it is very easy to find the matrix P in terms of the polynomial V . It is given 
by 



P 



dV„ 



all X, 



(5.9) 



where, as previously defined, indicates the term in V associated to the perfect 
matching p^. For our example, we obtain 
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(5.10) 



6. BFTs and Calabi-Yau's: Moduli Spaces 

A remarkable feature of BFTs, which is at the center of their special properties, is 
that perfect matchings extremely simplify the computation of their moduli space. The 
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Figure 20: The seven perfect matchings for the BFT in Figure [l8|. Edges in the perfect 
matchings are indicated in red. 



moduli spaces are automatically toric and perfect matchings are in one-to-one corre- 
spondence with fields in their gauged linear sigma model (GLSM) description. In this 
section we discuss this calculation, using the example in Figure |TB| to illustrate our 
ideas. Indeed, perfect matchings automatically satisfy F-term equations in the gauge 
theory. This property was originally identified in the context of dimer models on T 2 
in [|l| , and a detailed proof of how perfect matchings parametrize the moduli space of 
the corresponding theories was given in fr4lfl . We now briefly review the arguments in 
these papers, which extend without changes to general BFTs. 



6.1 F-Flatness and Perfect Matchings 



The map between chiral fields in the BFT and perfect matchings given in (2J) implies 
that F-term equations are trivially satisfied, as we now review. The vanishing of F- 
terms for fields associated to external legs is not imposed. 10 For any bifundamental 
field Xq associated to an internal edge, we have 



W = X P 1 (X l ) - X P 2 (X i ) + . . . , (6.1) 

where we have identified the only two terms in the superpotential containing Xq. Pi(Xi) 
and P2(Xj) are products of bifundamentals fields. The F-term equation for Xq takes 
the form 



d Xo W = « P 1 (X i ) = P 2 (X i ). (6.2) 

10 Since these fields appear in a single superpotential term, they would set to zero the product of 
fields they are coupled to. This special treatment of external legs is motivated by the connection to 
geometry, which we develop in this section. 
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This equation has a simple graphic representation as shown in Figure |2T]. After remov- 
ing X , the product of edges connected to node 1 needs to be equal to the product of 
edges connected to node 2. Using (|2.1| ), this becomes 

nn^=niK M ( 63 ) 

But this equation is automatically satisfied. Since nodes 1 and 2 are precisely separated 
by a single edge, every perfect matching that appears on the L.H.S. of ( |B~B|) also appears 
on its R.H.S. 

Pipy p,oy 

Figure 21: Graphic representation of the F-term equations in a BFT. 



6.2 The Master Space 

The first step in our discussion of the vacuum structure of BFTs is the concept of 



master space, which was introduced for arbitrary M = 1 field theories in [42). 11 The 
master space is defined as the space of solutions to F-term equations. Since D-terms 
are not imposed, we can regard the master space as the full moduli space of the gauge 
theory, including baryonic directions. 

Following the discussion in the previous section, the master space of a BFT is 
naturally parametrized in terms of perfect matchings. For a BFT, the master space is 
toric, i.e. it can be described by a GLSM. In GLSM language, F-term conditions can 
be translated to certain £7(1) charges of the perfect matchings, which are encoded in a 
charge matrix Qp defined as 

Qf = KerP. (6.4) 

The toric diagram of the master space is given by KerQp, which is indeed P. In other 
words, the matrix P connecting chiral fields in the quiver to perfect matchings gives 
the positions of points in the toric diagram of the master space! It is interesting to note 



l The concept of master space extends to SUSY theories in other dimensions. 
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that a few months before the general concept of master space was introduced in [421, 



the same object was constructed in the mathematics literature in f43| , from a different 
point of view, for the restricted case of bipartite graphs on a disk. In that work, the 
matrix P was referred to as the matching polytope. Our interpretation of the graph as 
defining a gauge theory makes the emergence of this geometry absolutely natural and 
allows its generalization to bipartite graphs on arbitrary Riemann surfaces. 

As we have just said, the toric diagram for the master space is given by the P 
matrix. In order to obtain a better idea of this geometry, it is useful to consider the 
row-reduced version of P which, for the example at hand, becomes 
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(6.5) 



We conclude that the master space is a 5-complex dimensional toric geometry with a 
toric diagram consisting of seven different points. From now on, every time we mention 
the dimension of a Calabi-Yau manifold, we refer to its complex dimension. Further- 
more, the entries in every column of G mast add up to 1, implying the master space is 
Calabi-Yau. In fact, the Calabi-Yau property will be exhibited by the master spaces 
of all models considered in this paper. Since the toric diagram lives on a hyperplane 
at distance 1 from the origin, we can project it down to 4 dimensions by, for example, 
considering only four of the rows in ( |6.5| ). A convenient way of visualizing this 4d toric 
diagram is by considering different 3d projections, as shown in Figure [221. Different 



points in the 5d toric diagram might be projected down to the same point in 3d. Such 
points are indicated in red in Figure |2|. 

The interior of the toric diagram of the master space, i.e. of the matching polytope, 
provides a graphical representation of the corresponding cell in the positive Grassman- 
nian. The BFT interpretation of the lower dimensional sub-cells on its boundary will 
be discussed in Section [K| 

Taking the kernel of P, or equivalently of G mast , we obtain the charge matrix that 
implements the F-terms 
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(6.6) 
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(a) 



(b) 



Figure 22: Two projections of the toric diagram for the master space, corresponding to 
( |6,5[ ). Points descending from multiple ones in 5d are shown in red. The projections cor- 
respond to keeping the following combinations of rows: a) {G rri ast,i-,G rri ast,2-,G mas tfi) and b) 

{G m ast,2i G m ast,3i G m ast,A) • 



6.3 The Mesonic Moduli Space 

The mesonic moduli space, is another natural geometry associated to any BFT. For 
shortness, we will refer to it simply as the moduli space from now on. The moduli space 
of any gauge theory is the vacuum space of solutions of both vanishing F and D-terms. 
It is thus a projection of the master space onto the subspace of vanishing D-terms. 
There is a D-term contribution for each gauge group in the BFT i.e., by means of 



the dictionary introduced in Section |2T|, for every internal face in the bipartite graph. 
It is convenient to define the charge matrix A of the BFT, as the matrix encoding 
how every chiral field transforms under the gauge symmetries. 12 The matrix A is an 
n fields x rigauge = n edges x n int . faces matrix in which rows correspond to chiral fields 
and columns correspond to gauge groups. For the row associated to X^-, the non-zero 
entries are a 1 for the i th column and a —1 for the j th column. All entries are zero in 
rows associated to adjoint fields X^. D-terms can then be encoded in a charge matrix 
Qd giving the charge of perfect matchings under the gauge groups. This means that 
Qd is defined such that 



P-Q 



A. 



(6.7) 



It is clear that (3/7) does not determine Qd uniquely. Any solution to this equation is 
equivalent for the purpose of determining the moduli space. 



2 For the purpose of this paper, it is sufficient to proceed as if every gauge group was U(l). 
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For our example, there is a single gauge group associated to face 1. Under it, X i3 
and X 15 have charge 1, X 2 \ and X 41 have charge —1, and all other fields are neutral. 
It is straightforward to verify that the following Q D does the right job 



Q 



D 



Pi P2 P3 Pi Pb P% Pi 

0-11-11 



(6i 



The next step in the determination of the moduli space is to concatenate Qf and 
Qd into a single charge matrix Q 

« = (&)■ (M) 

The toric diagram of the moduli space is thus encoded in a matrix G such that 

G = KerQ. (6.10) 
Let us consider our example. From ( |6.6|) and ( |6.8|) , we obtain 
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(6.11) 



The moduli space is a 4d toric manifold. As for the master space, the entries in 
every column add up to 1, implying the moduli space is also a Calabi-Yau manifold. 
This will also be the case for all the examples considered in the paper which, together 
with our previous observation regarding the master space, leads us to conjecture that 



Conjecture: 

The master and moduli spaces of every BFT are toric Calabi-Yau manifolds. 



We expect the existence of a simple proof of this statement based on the combinatorics 
of Qf and Qd- 

We observe an interesting phenomenon: there can be non-trivial multiplicities of 
perfect matchings associated to the same point in the toric diagram. In particular, 
we see that the point (—1, 1,0, 1) corresponds to both pi and p 2 - Such multiplicities 
are generic in BFTs. For example, for the case of D3-branes probing toric CY 3-folds, 
trying to understand them has been an important factor leading to the correspondence 
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between the associated quivers and dimer models [44|. The role of multiplicities in the 
generalized context of BFTs is certainly an interesting question that deserves further 
investigation. 

Multiplicities can arise for both internal and external points of a toric diagram. We 



would like to note that the toric diagram in Figure |23] contains corners with multiplicity 
different from one. In the specific case of gauge theories on D3-branes probing toric 
CY 3- folds (i.e. BFTs on T 2 with no boundaries), this feature has been identified as 
an indication of an inconsistency. In fact, for this class of theories, this behavior is 
directly connected to the existence of self- intersecting zig-zag paths |22]. This is clearly 
not the case here, as we can verify by explicit determination of the zig-zag paths. At 
this time, we are not aware of any pathology signaled by this behavior. 13 In any case, 
we are confident that BFTs without this feature exist. Such models can be analyzed 
with exactly the same methods we have applied here. 

Since the toric diagram lives on a hyperplane at distance 1 from the origin, it can 
be projected down to three dimensions by, for example, considering any three of the 
rows in (6.11). Figure shows the toric diagram for this model. This is indeed a 
well-known geometry, the real cone over the 7-dimensional Sasaki-Einstein manifold 
Q 1,1,1 . This CY 4- fold has been extensively investigated in connection to other types 
of gauge theories associated to M2-branes J|5|, [|(| |47|, |48| . 



Figure 23: Toric diagram for the CY 4-fold that is the moduli space of the 4-leg, 1-loop 
model. This CY 4-fold is the real cone over Q ' 4. We indicate the non-trivial perfect 
matching multiplicity of the (—1, 1,0, 1) point in the toric diagram with a number. 



The moduli space of a gauge theory is invariant under Seiberg duality. More 
abstractly, in graph theoretic language, this means that the moduli space of a BFT 

13 This diagnostic quite probably does not apply to other classes of theories beyond those on T 2 . 
For instance, there are known examples of gauge theories in 2+1 dimensions whose moduli space have 
toric diagrams with corner multiplicities and that do not have any known problem G7J. 
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provides a natural geometry associated to a bipartite graph on a Riemann surface that, 
by construction, is invariant under square moves. In Section §, we will discuss in detail 
the implications and applications of this fact. 



7. Additional Examples: Increasing Boundaries and Genus 

BFTs associated to graphs without boundaries have been extensively studied in the 
literature. Tilings of T 2 describe the gauge theories on D3-branes over toric CY 3-folds 
JTJ, and tilings on higher genus Riemann surfaces arise when acting on them with 
the untwisting map as discussed on Section |] ||. For this reason, we emphasize in 
this section the novel case of BFTs with boundaries, which are also the ones that are 
relevant for scattering amplitudes. We start discussing a model on the disk and soon 
move to theories that have never been studied before: models with multiple boundaries 
and higher genus. It is natural to expect such configurations to be relevant for leading 
singularities beyond the planar limit. Further studies of non-planar graphs will appear 



in 



We will put special emphasis in the geometry of the corresponding master and 
moduli spaces. They can be determined in terms of perfect matchings following the 
general procedure introduced in Section 0. 



7.1 Another Example on the Disk: The Hexagon-Square Model 



Let us consider the 2-loop graph, shown in Figure |24], corresponding to the scattering of 
3 negative helicity and 3 positive helicity gluons. For pedagogical reasons, we present 
the full details of the calculation of its master and moduli spaces in Appendix [X]. The 
treatment of other examples in the paper will be briefer and will only emphasize the 
main results. The master Kasteleyn matrix for this model is 
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(7.1) 



The theory has 25 perfect matchings. The master space is an 8d toric CY. Its toric 
diagram is given by the matrix 
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(a) (b) 

Figure 24: a) Bipartite graph for the hexagon-square model. It contains two internal and 
six external faces, b) The six zig-zag paths for this model. We see that it does not have any 
self-intersecting zig-zag path. 
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The moduli space is a 6d toric CY. The 25 perfect matchings give rise to 18 different 
points in its toric diagram, whose positions are captured by the following matrix 
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(7.3) 



where we have introduced a compact notation for the G matrix, in which the last 
row indicates the multiplicity of perfect matchings for each point in the toric diagram. 
As already done in Section ^ a useful way of visualizing this 6d toric diagram is by 



considering different 3d projections, as shown in Figure 25. 
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(a) (b) 

Figure 25: Two projections of the toric diagram corresponding to ( |7.3|) , Points descending 
from multiple ones in 6d are shown in red. The numbers indicate the non-trivial multiplicity 
of perfect matchings. The projections correspond to keeping the following combinations of 
rows: a) (Gi - G 2 + G 3 , G 3 , G 4 + G 6 ) and b) (G U G 2 , G 3 ). 



7.2 Two Boundaries: the Cylinder 



We now study a model with more than one boundary. Let us consider the example in 
Figure It has 6 external legs distributed on 2 boundaries, and 3-loops. 




The master Kasteleyn matrix is 
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This model has 44 perfect matchings. The master space is a lOd CY. The moduli 
space is a 7d CY. The 44 perfect matchings organize into 28 distinct points in the toric 
diagram, which are given by the following matrix 
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Figure |27] shows two possible 3d projections of this toric diagram. 




Figure 27: Two projections of the toric diagram corresponding to ( [7.5[ ). Points descending 
from multiple ones in 7d are shown in red. The numbers indicate the non-trivial multiplicity 
of perfect matchings. The projections correspond to keeping the following combinations of 
rows: a) (G 1 ,G 2 ,G 3 ) and b) (G 2 - G 3 ,G 4 ,G 5 ). 
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7.3 One Boundary on T 2 



Let us move to higher genus and consider a model on a 2-torus, with 4 external legs 
terminating on a single boundary. From a scattering amplitude perspective, we can 
regard this diagram as a non-planar, 5-loop contribution to the scattering of 2 negative 
helicity and 2 positive helicity gluons. We refer to this theory as model 1 and we show 
the corresponding graph in Figure . 
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Figure 28: a) Bipartite graph for model 1. It lives on a 2-torus and has four external nodes 
on a boundary. It contains five internal and four external faces, b) The four zig-zag paths 
for this model. 



The master Kasteleyn matrix for this model is 
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The theory has 48 perfect matchings and the master space is an lid CY. The 
perfect matchings give rise to 22 different points in the toric diagram of the moduli 
space, with positions summarized by the following matrix 



35 



/ 1 1 1 



10 

10 



111 



1-10 
10 1 



110 1 




1 
1 



1 -1 1 -1 -1 1 
01-1-11 000 2-2 



-1 -1 1 -1 -1 1 2 \ 

1-10 1-10 0-11 
10 110 10 10 
110 110 10 

1 1 1 -1 -1 -1 
-1 -1 -1 1 1 1 



\ 9 4 4 3 3 33322 1 1 1 1 1 1 1 11111/ 



We see that the moduli space is a 6d toric CY. In Figure [29|, we show two projections 
of the toric diagram down to three dimensions. They have been chosen in order to 
minimize the overlap of distinct points after the projections. In both cases, the 22 
points of original 6d toric diagram are mapped to 19 points. 




Figure 29: Two projections of the toric diagram corresponding to (|7.7|). We show in red the 
points descending from multiple ones in 6d. The numbers indicate the non-trivial multiplicity 
of perfect matchings. The projections correspond to keeping the following combinations of 
rows: a) (Gi,G 2 ,G 6 ) and b) (G U G 2 ,G 5 - G 3 ). 



8. Square Moves and Geometry, or Seiberg Duality and Moduli 
Spaces 

Let us investigate the effect of square moves on BFT theories. As we have explained in 
Section j|, they correspond to Seiberg dualities on certain gauge groups of the BFTs. 
The moduli space of the theories is, by construction, invariant under Seiberg duality. 14 

14 As shown in this is not the case for the master space which, in the case of plabic graphs, is 
the toric geometry associated to the matching polytope Q. 
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As a result, the moduli space is an ideally suited object for identifying theories con- 
nected by square moves. This problem becomes rather non-trivial for large graphs, 
multiple square moves, multiple boundaries and/or higher genus Riemann surfaces. 



8.1 The Dual of the Hexagon-Square Model 

Let us consider the model shown in Figure RO, which is obtained from the hexagon- 



square model discussed in Section [T]] by Seiberg dualizing the gauge group associated 
to face 2. Four 2-valent nodes, i.e. mass terms in the BFT, are generated by the duality. 
We have only integrated out the massive fields associated to two of them, in order to 
preserve the external legs connected to nodes 6 and 13. 



8 



13 




14 

(a) 



(b) 



Figure 30: a) Bipartite graph obtained by Seiberg dualizing the gauge group associated to 
face 2 of the hexagon-square model. It contains two internal and six external faces, b) The 
six zig-zag paths for this model. 



The master Kasteleyn matrix is 
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The theory has 22 perfect matchings and the master space is an 8d toric CY. The 
moduli space is a 6d CY, with toric diagram given by 
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This moduli space is the same as the one for the original hexagon- square model. It is 
indeed possible to find an SX(6, Z) transformation that takes ( |S.2|) into ( [7.3|) . Instead 
of giving the explicit transformation, we show two projections of the toric diagram in 



Figure 31, which are identical to those shown in Figure 25 for the Seiberg dual theory. 





(b) 



(a) 

Figure 31: Two projections of the toric diagram corresponding to ( |3,2D . Points descending 
from multiple ones in 6d are shown in red. The numbers indicate the non-trivial multiplicity 
of perfect matchings. The projections correspond to keeping the following combinations of 
rows: a) (G a - G 4 , G 5 , G 6 ) and b) (G 2 - G 3 , G 4 , G 6 ). 



Comparing ( \l.3\) and Q3.2|) , we see that the original hexagon- square model and its 
Seiberg dual differ in the multiplicity of perfect matchings associated to each point in 
the toric diagram of the moduli space. This is a generic feature of dual theories that will 
be also encountered in the examples that follow. Different multiplicities are a manifes- 
tation, in the context of toric geometry, of the action of cluster transformations. They 
relate the partition functions for perfect matching of dual models and, in particular, 
produce the perfect matching multiplicities for any of the two theories in terms of those 
for the other one. The central role of cluster transformations, which leave the boundary 



measurement invariant [26[] , in the study of leading singularities has been investigated 
in M. A more intuitive understanding of the role of cluster transformations in the BFT 



will be presented in future work f52fl . 



8.2 Non-Planar Duals 

Let us now consider a theory, that we call model 2, obtained from model 1 in Section 
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[7.3| by Seiberg dualizing face 4. The resulting graph is shown in Figure 
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(a) (b) 



Figure 32: a) Bipartite graph for model 2, obtained from model 1 by Seiberg dualizing face 
4. It lives on a 2-torus and has four external nodes on a boundary. It contains five internal 
and four external faces, b) The four zig-zag paths for this model. 

The master Kasteleyn matrix is 
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The theory has 53 perfect matchings, which turn into 22 distinct points in the toric 
diagram of the moduli space. We already see that this number agrees with its Seiberg 
dual. The moduli space is a 6d toric CY whose toric diagram is given by 
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This moduli space is identical to the one for model 1. As in the previous exam- 
ple, instead of providing the explicit SL(6, Z) transformation connecting the two toric 
diagrams, we present some 3d projections in Figure |33], which match those in Figure |29. 




Figure 33: Two projections of the toric diagram of model 2, given by (|3.4|). Points descending 
from multiple ones in 6d are shown in red. The numbers indicate the non-trivial multiplicity 
of perfect matchings. The projections correspond to keeping the following combinations of 
rows: a) (Gi, G2, Gq) and b) {G\, G2, G5 + G3). Models 1 and 2 have the same moduli space. 

Despite their simplicity, the examples considered in this section show how powerful 
the concept of moduli space is for identifying models connected by Seiberg duality, 
i.e. configurations related by square moves in the graph. The moduli space serves as 
a practical and sharp diagnostic even for large graphs, complicated topologies and/or 
theories related by a chain of multiple Seiberg dualities. 



9. Loop Reduction and Calabi-Yau Geometry 

In this section we study some nice behavior exhibited by the moduli spaces associated 
to multi-loop diagrams. For concreteness, we consider the scattering of 2 negative and 2 



positive helicity gluons and focus on the multi-loop ladder diagrams given in Figure |34|, 
which generalize the model studied in Section [| 15 

For BFTs on a disk, the dimension of the master space is equal to the total number 
of faces of the graph []43" |. To determine the moduli space, we further impose the D-term 



equations associated to internal faces, so the dimension of the moduli space is equal to 
the number of external faces. Applying this general discussion to the class of models 



given by Figure 34, we conclude that while the master space of the n-loop theory is 



15 



This example was also independently considered by the authors of 
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Figure 34: The n-loop ladder diagram for the scattering of 2 negative and 2 positive helicity 
particles generalizing the 1-loop model discussed in Section 0. 



(n + 4)-dimensional, the moduli space is a CY 4-fold for every n. We will soon see that 
the agreement between moduli spaces goes beyond just the number of dimensions. 

Two Loops 



Let us begin with the 2-loop diagram shown in Figure |35[ The master Kasteleyn matrix 
for this model is 
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/ 



The model has 10 perfect matchings and the master space is, as anticipated, a toric 





(a) (b) 

Figure 35: a) Bipartite graph for the 4-leg, 2-loop model. It contains two internal and four 
external faces, b) The four zig-zag paths for this model. 
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6d CY. This theory has two gauge groups, associated to faces 1 and 2. Performing the 
further quotient by these symmetries we obtain the moduli space, which is a toric 4d 
CY. Its toric diagram is given by the matrix 



G 



( 1 -1 \ 
1 10 1 
10 1 10 
-1 1 
V 3 2 2 111 j 



(9.2) 



The associated toric diagram is presented in Figure [36| Interestingly, the moduli 
space is exactly the same as for the 1-loop model considered in Section ||, although 
with different perfect matching multiplicities (a, b, c, d, e, /) = (2, 2, 1, 1, 1, 3). 
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Figure 36: General toric diagram for the moduli space of the 4-leg, 2, 3 and 4-loop models. 
The moduli space is in the three cases the same as the one for the 1-loop model, i.e. the real 
cone over Q ' ,1 . Letters indicate the non-trivial perfect matching multiplicities of points in 
the toric diagram, which depend on the number of loops. 



Three Loops 



Let us now quickly analyze the 3-loop model, given by the graph in Figure [J7]. The 
master Kasteleyn matrix is 
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(a) (b) 

Figure 37: a) Bipartite graph for the 4-leg, 3-loop model. It contains three internal and 
four external faces, b) The four zig-zag paths for this model. 



The model has 15 perfect matchings and a master space that is a CY 7-fold. After 
imposing the D-terms associated to the three gauge groups, we see that the moduli 
space is a CY 4-fold with toric diagram given by 
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We conclude that the moduli space of the 3-loop graph is identical to the one for 1 
and 2-loops, shown in Figure [36], with multiplicities (a, b, c, d, e, /) = (5, 2, 1, 3, 1, 3). 

Four Loops 



The 4-loop model corresponds to the graph in Figure |38|. The master Kasteleyn matrix 
is 
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The model has 23 perfect matchings and a master space that is a CY 8-fold. The 
moduli space is a CY 4-fold with toric diagram given by 
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Figure 38: a) Bipartite graph for the 4-leg, 4-loop model. It contains four internal and four 
external faces, b) The four zig-zag paths for this model. 
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which, once again, precisely agrees with the moduli spaces for the lower loop models, but 
with multiplicties (a, b, c, d, e, f) = (8, 3, 1, 5, 1, 5). 

It is natural to expect that the moduli space remains the same, up to perfect matching 
multiplicities, for arbitrary number of loops. Below we show that this is indeed the case 
by proving that different loops are connected by Seiberg duality. We conjecture that this 
behavior is a geometric manifestation of the fact that for a given set of scattered particles 
the number of leading singularities is finite and determines the scattering amplitude to an 
arbitrary number of loops. It would be interesting to investigate this phenomenon for other 
multi-loop diagrams. 



9.1 Loop Reduction and Seiberg Duality 

The fact that BFTs associated to the ladder diagrams with an arbitrary number of loops 
share the same moduli space suggests that they are connected by Seiberg duality. Since a 
different number of loops maps to a different number of gauge groups in the BFT, Seiberg 
duality clearly needs to be supplemented with some additional dynamics, as we now explain. 
The n-loop diagram is connected to the (n — 2)-loop one by the sequence of steps shown 



in Figure 39. In terms of gauge theory dynamics, these steps have the following meaning: 



1) We perform a Seiberg duality on some of the internal faces that are not at the endpoints 
of the ladder. This transformation is implemented by a square move in the graph and 
generates four 2-valent nodes. 
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Figure 39: An n-loop diagram can be turned into an (n — 2)-loop one by a sequence of steps 
that in the BFT correspond to: 1) Seiberg duality, 2) Integrating out massive fields, and 3) 
Confinement of Nf = N c gauge groups and formation of mesons. 



2) The 2-valent nodes generated in the previous step correspond to mass terms in the 
superpotential. We integrate out the massive fields, which maps to condensation of 
nodes in the graph. When doing so, the number of sides of each of the two faces 
adjacent to the dualized one is reduced to two. 

3) Internal faces with two sides correspond to SU(N C ) gauge groups with Nf = N c . At 
low energies, such gauge groups confine and their dynamics is expressed in terms of 
gauge invariant (under the confined gauge group) mesons and baryons. The graphic 
implementation of the formation of mesons corresponds to combining the two edges on 
the boundary of these faces into a single one. This process makes the faces disappear, 
in agreement with confinement. 

Iterating this process we can show that all diagrams with an even number of loops give rise 
to dual gauge theories. Similarly, theories with an odd number of loops are also dual. 



10. The Boundary Operator as Higgsing 



In Section 4.4, we reviewed the correspondence between bipartite graphs and cells in the 



positive Grassmannian. In Section p.2| we explained, in terms of the master space, that a cell 
takes the form of a convex polytope and its boundary is a collection of lower dimensional cells. 
Being at a boundary of a cell corresponds to setting some of the entries of the corresponding 
matrix C to zero. The larger the number of vanishing entries is, the lower dimensional the 
corresponding boundary cell is. 
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Setting an entry in C to zero corresponds to eliminating the connectivity between the 
associated external nodes. This is achieved by removing an internal edge in the graph, 



disrupting oriented paths between the nodes. Figure 40 shows an example of this process. 
Removing the edge shown in red results in setting C24, C26, C27, C37 and C57 to zero in (4.3). 



The discussion in Section 3.3 implies that the boundary operator maps to higgsing in the 
BFT. 




Figure 40: An example of higgsing. An edge is removed from the graph, resulting in the 
disappearance of some oriented paths, associated to a perfect orientation, connecting external 
nodes. 



Turning on a non-zero vev for a bifundamentals field Xq determines an energy scale 
{Xq) = A. Removing the edge associated to Xq from the graph corresponds to considering 
energies much smaller than A. The surviving graph accurately captures the low energy 
physics, such as the moduli space of the theory, provided the vevs involved are much smaller 
than A. 

The BFT not only captures the combinatorics of the boundary but also describes the 
continuous approach to boundary facets. In fact, identifying expectation values of bifunda- 
mental fields with the inverse of edge weights entering (}4.2| ), we obtain perfect agreement 
between the higgsing and Grassmannian pictures. 16 As some expectation value is increased 
certain entries in C get suppressed, eventually vanishing once the vev is sent to infinity. 

We have discussed how the boundary operator is linked to a simple local operation on 
the graph: edge removal. Given the map between zig-zag paths and Deligne permutations 



explained in Section 2.3, it is straightforward to see that the boundary operator acts by 



16 



The inversion of edge weights in this correspondence is not surprising. The structure of F-terms 
in BFTs is such that there is a trivial x e x~ x symmetry. 
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flipping Deligne permutations §, according to 



fd(a) = b 
fd(c) = d 







fd(a) 
fdic) 



d 
b 



(10.1) 



where a, b, c and d are the endpoints of two intersecting zig-zag paths. Figure fil] shows how 
the permutation flip results from edge removal. Notice that while permutation flip is a global 
operation in the graph, it is equivalent to removing edges, which is local. 




fM = b 



Figure 41: Removing an internal edge gives rise to a recombination of zig-zag paths that 
results in a flip of Deligne permutations. 



10.1 Consistent Higgsing and Untwisting 



We should only consider higgsings that produce consistent graphs, i.e. graphs without self- 
intersecting zig-zag paths. When an edge is removed, two zig-zags are split at some intermedi- 



ate points and then recombined as in Figure 41. This implies that removing an edge generates 
self-intersections only when the two zig-zags involved originally have multiple intersections. 



This situation is sketched in Figure 42 




Higgsed 



Figure 42: Removing an edge between two zig-zag paths with multiple crossings results in 
a zig-zag with multiple intersections. 



Zig-zags are not manifest in the graph G, and keeping track of them or recomputing them 
after each higgsing in order to check consistency is rather tedious. This is particularly hard 
when multiple non-zero vevs are involved. It is then useful to consider the untwisted graph 
G. Zig-zags of G become boundaries of faces (both internal and external) in G and explicit 
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in the graph even after higgsing. Recombination of zig-zags maps to the recombination of 
faces and we can efficiently follow them through the process of removing edges. 

According to our previous discussion, in order to preserve consistency, edges between 
zig-zags that intersect more than once cannot be removed. It is then straightforward to 
identify inconsistent higgsings using G: we simply cannot delete edges sitting between faces 
with multiple intersections. 

Let us illustrate these concepts with an explicit example. Consider Figure which 
corresponds to a leading singularity in the scattering of 4 negative and 4 positive helicity 
gluons at 4-loops. We have labeled edges to facilitate their identification after untwisting. 




(a) (b) 

Figure 43: a) A 4- loop diagram associated to the scattering of 4 negative and 4 positive 
helicity gluons. We have labeled edges in blue, b) The eight zig-zag paths of this model. 

There are eight zig-zag paths: two of length 4, two of length 6, two of length 7 and two 
of length 8. They are given by the following collections of edges: 



Length 4 1 : (1, 4, 10, 14) Length 7 3 : (1, 5, 11, 17, 18, 19, 20) 

2 : (2, 8, 13, 20) 4 : (2, 7, 12, 17, 16, 15, 14) 

Length 6 5 : (3, 10, 15, 21, 24, 25) Length 8 7 : (3, 4, 5, 6, 12, 18, 22, 25) 

6 : (9, 13, 19, 22, 24, 23) 8 : (9, 8, 7, 6, 11, 16, 21, 23) 



(10.2) 



The untwisted graph G is shown in Figure [W|. 18 The eight zig-zag paths of the original 
graph turn into eight external faces via untwisting. G has three boundaries and lives on a 
2-torus. 



17 It is important to emphasize that we do not need to require G to be consistent. In fact, it has 
self-intersecting zig-zags if the original theory has adjoint fields. 

18 We are extremely grateful to Rak-Kyeong Seong for sharing his expertise in untwisting complicated 
graphs and for verifying this example. 
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Figure 44: The graph G obtained by untwisting the zig-zag paths in Figure Face labels 
are those of the corresponding zig-zag paths in the original graph. G contains eight external 
faces, has three boundaries (represented as blue circles) and lives on a 2-torus. 



As explained above, G can be used to identify, by simple inspection, all edges that cannot 



be removed without spoiling consistency. They are indicated in color in Figure 44. In green, 
we show external legs, which correspond to the scattered particles and hence are preserved. 
In red, we show the edges that cannot be removed because they sit between faces with more 
than one intersection. In particular we have edges 5 and 18 between faces 3 and 7, and edges 
7 and 16 between faces 4 and 8. While this approach might seem too elaborate for dealing 
with the removal of single edges, it becomes particularly useful for systematically determining 
consistent removals of multiple edges. 

10.2 Higgsing and Geometry 

Let us discuss the effect of higgsing on the geometry of the master and moduli spaces. The 
correspondence between perfect matchings which should be regarded as GLSM fields, and 
a chiral field X, in the quiver is given by (2.1), which we reproduce here for convenience 



X l = ]Jp^. (10.3) 

This means that, in order for Xi to acquire a non-zero vev, all perfect matchings with Pj„ = 1, 
i.e. all the ones that contain the corresponding edge in the graph, need to get a non-zero 
vev and are eliminated. This removal of perfect matchings results in the disappearance of 
some points in the toric diagrams of the master and moduli spaces, and in a reduction in the 
multiplicity of others. 

While higgsing removes points from the toric diagrams of the master and moduli spaces, 
their qualitative behavior is different. Higgsing decreases the number of internal faces in G 
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by one, which implies that the dimensionality of the master space is reduced by the same 
amount. The toric diagram of the master space of the higgsed theory is, as expected, a lower 
dimensional sub-cell on the boundary of the original one. On the other hand, the number of 
D-term equations is also reduced by one when higgsings, which implies that the dimension of 
the moduli space remains constant. 



11. Conclusions and Outlook 

In this paper, we have introduced and started the study of Bipartite Field Theories, a general 
class of 4d, Af = 1 quiver gauge theories defined by a partite graph on a Riemann surface. We 
explained the map between the field theory dynamics and graph modifications, the emergence 
of toric CY manifolds as the master and moduli spaces of the gauge theory, the connection 
between Seiberg duality and square moves and interpretation of the boundary operator on 
cells in the positive Grassmannian as higgsing. While our discussion has been completely 
general, most of our explicit examples have boundaries, i.e. they were relevant for scattering 
amplitudes. 

We developed specific tools, in the form of generalized Kasteleyn matrix techniques, for 
the study of general BFTs. We extensively demonstrated in examples how these methods 
allow the explicit determination of master and moduli spaces, even for complicated graphs. 
Finding the corresponding CY manifolds is reduced to computing determinants of rather 
simple matrices. 

We have studied models that go beyond the planar limit, by increasing the number of 
boundaries and the genus of the underlying Riemann surface. From a BFT viewpoint, they 
are not more involved than planar graphs. As discussed in Section || it is in the context 
of non-planar graphs that the moduli space of the BFT more strikingly shows its power for 
identifying theories related by square moves, since they cannot be understood in terms of 
permutations as in the planar case. 

Our work suggests various interesting directions for future investigation, some of which 
are summarized below: 



One of the main lessons of our work is the existence of a deep connection between 
Calabi-Yau manifolds, which appear in the form of master and moduli spaces of BFTs, 
and leading singularities. As explained in Sections 3.2 and 1C, the master space is the 
natural object for describing cells in the Grassmannian and their boundaries. But, to 
us, the moduli space appears to be a more fundamental object, since it is invariant under 
all equivalence transformations of the graph. It is natural to expect these geometries to 
play central role in the study of leading singularities, yet to be unveiled in full generality. 
In the restricted case of graphs on a disk, the toric diagram of the moduli space exhibits 
striking similarities with the matroid polytope of jSJ. The connection between these 
two objects will be explained in |p0|. 



50 



For graphs without boundaries on T 2 , there is an intimate connection between zig-zag 
paths and the geometry of the moduli space of the associated BFT, which is a CY 3- 
fold. Zig-zags are in one-to-one correspondence with (p, g)-legs normal to the faces of the 
toric diagram. It would be interesting to understand whether, and if so how, zig-zags on 
bipartite graphs on generic Riemann surfaces with boundaries are related to the moduli 
spaces of the corresponding BFTs. We expect that understanding this connection, in 
conjunction with the previous point, will lead to a map between invariants, global 
properties of zig-zags and CY manifolds. 

A string theory embedding is known for certain classes of BFTs. This is the case for 
those associated to graphs without boundaries on T 2 , which arise on the world volume 
of D3-branes probing toric CY 3-folds. In these cases, the graph can be interpreted 
as a physical web of NS5-branes from which D5-branes are suspended, connected by 
two T-dualities to the original configuration of D3-branes. In addition, the graph on a 
(typically) higher genus Riemann surface that is obtained from it by the untwisting map 
describes a configuration of D6-branes on the mirror manifold. It would be extremely 
interesting to investigate whether additional sub-classes of BFTs admit a string theory 
realization. 

In Section ||, we discussed how different BFTs can give rise to the same CY manifolds 
as their moduli spaces. We can regard such equivalence as a generalized version of 



Toric Duality |53|, |54], [55|, ^6|. In Section g, we presented explicit examples showing 
that even graphs with a different number of loops can lead to the same CY. There 
is a priori no obvious reason preventing graphs with different numbers of boundaries 
and/or genus from producing the same CY. 19 It would be interesting to search for 
explicit examples realizing this behavior and, if they indeed exist, to understand what 
its physical interpretation is. 

The Grassmannian G(k, n) also arises as the moduli space of k U(n) vortices. Given 
the connection between cells in the Grassmannian and bipartite graphs and the fact 
that any such graph can be reduced to one only involving 2 and 3-valent nodes, it would 
be interesting to explore whether some natural decomposition of the moduli space, in 
which U(2) and U(3) vortices are basic building blocks, exist. 

We foresee BFTs will provide useful insights and guidance and will fuel future devel- 
opments in the study of systems associated to bipartite graphs, in particular in the area of 
leading singularities in scattering amplitudes. 



19 In fact we know that a similar behavior is possible for the master space: specular duals, whose 
underlying graphs are related by untwisting zig-zag paths, share the same CY as the master space 
31 . 
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A. Master and Moduli Spaces of the Hexagon- Square Model 



The bipartite graph for this model is shown in Figure 24. The master Kasteleyn matrix for 
this model is 

/ 





8 


9 


10 


11 


12 


13 


14 


1 


X31 








Xis 


A~83 








2 


Xu 


X42 


X21 














3 





X25 


Xq2 








X56 





4 








Xie 


Xji 








X&7 


5 


X43 




















6 





X54 

















7 











X87 












(A.l) 



from which we determine 

V = —^14^16^18^25 — ^16^18^42^43-^56 + ^14^16^18^54^56 _ -^18^21^25^43^67 

— A18X21X43X54A56A67 + Xx^X^X/^X^X^ — X14X18A54A62A67 — X21X25X31 X71 

+X21X31A54A56 A71 + X3iXl2X62^71 + ^21X25X43X71X33 — X2 1X43X54X56X71X33 

— X42X43X62X71X83 + X14X54X62X71X83 + X16X31X42X56X87 + X21X25X31X67X87 

— X2iX3iX 5 4X 56 X 6 7X87 — X 3 iX42Xg2X 6 7X 8 7 + X 14 XigX25X83X87 — Xi 6 X42X43X 56 X83X87 

+X14X16X54X56X83X87 — X21X25X43X67X83X87 + X21X43X54X56X67X83X87 

+X42X43X62X67X83X87 — X14X54X62X67X83X87. (A. 2) 



I.e. the model has 25 perfect matchings. The P matrix becomes: 
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Pi 


P2 


P3 


Pi 


P5 


P6 


P7 


P8 


P9 


P10 


Pll 


P12 


f>13 


P14 


P15 


P16 Pl7 PIS Pl9 P20 P21 P22 P23 P24 P21 
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1 


1 


1 


1 


1 


1 
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Xie 
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1 


1 
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1 
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1 
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1 


1 


1 
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1 
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1 
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1 
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1 


1 


1 


1 
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1 
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1 


1 


1 


1 
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1 





1 
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1 





Xd2 
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1 
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1 

















1 


1 
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1 
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1 
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1 
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1 





1 














1 


1 








1 


1 








1 
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V -^87 








1 


1 





1 








1 


1 











1 





1 





1 





1 





1 





1 






(A.3) 



We can get a better idea of the toric diagram of the master space by considering the row- 
reduced version of this matrix 



/ pi 


P2 


P3 


P4 


P5 


P6 


P7 


P8 


P9 


P10 


Pll 


P12 


P13 


P14 


P15 


P16 


P17 


P18 


P19 


P20 


P21 


P22 


P23 


P24 P25 


1 
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-1 











-1 
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1 
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—2 
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-1 




















1 























1 





1 





1 





1 





1 
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-1 -1 
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1 1 
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1 





1 


1 


1 


1 


1 














1 


1 


1 1 



(A.4) 



The master space is hence an 8d toric CY. The charge matrix encoding F-terms is Qf = 
KerP, which becomes 
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Qf = 



I pi 


P2 


P3 
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P6 


P7 


P8 


P9 


PlO 
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/ 



(A.5) 



This theory has two gauge groups. The matrix associated to D-terms can be chosen to 

be: 



/ PI P2 P3 Pi P5 P6 P7 P8 P9 PlO Pll P12 P13 P14 P15 P16 P17 P18 P19 P20 P21 P22 P23 P24 P25 \ 

Q D = 000000100 1 1 -1 -1 -1 . 
\00000000000000001000 -1 0000/ 

(A.6) 



Combining Qp and Qr> into the total charge matrix Q, we obtain the matrix G = Ker Q 
defining the toric diagram of the moduli space 



/ Pl P2 P3 P4 P5 P6 P7 P8 P9 PlO Pll Pl2 Pl3 P14 Pl5 P16 Pl7 Pl8 Pl9 P20 P21 P22 P23 P24 P25 \ 

1 -1 -1 -1 -1 -1 -2 -1 -1 -2 -1 -1 -1 1 
0011121 -1 00000101010101010 
G= 0011011010000110011001100 
1010111000110111111110000 
-1 -1 -1 -1 1 1 1 1 
\0000 -1 -1 -1 111 000000000000000/ 

(A.7) 

In (ffl]), we presented a reduced version of this matrix in which we only show the different 
column vectors and indicate their multiplicity. 
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